MASSLESS NEUTRINO OSCILLATIONS 



F. Benatti 

Dipartimento di Fisica Teorica, Universita di Trieste 
Strada Costiera 11, 34014 Trieste, Italy- 
arid 

Istituto Nazionale di Fisica Nucleare, Sezione di Trieste 



R. Floreanini 

Istituto Nazionale di Fisica Nucleare, Sezione di Trieste 
Dipartimento di Fisica Teorica, Universita di Trieste 
Strada Costiera 11, 34014 Trieste, Italy 



Abstract 

Quantum dynamical semigroups provide a general framework for study- 
ing the evolution of open systems. Neutrino propagation both in vac- 
uum and in matter can be analyzed using these techniques: they allow 
a consistent treatment of non-standard, dissipative effects that can al- 
ter the pattern of neutrino oscillations. In particular, initially massless 
neutrinos can give rise to a nonvanishing flavour transition probability, 
involving in addition the Majorana CP- violating mixing phase. 
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1. INTRODUCTION 



Elementary particle systems are usually treated as isolated quantum systems: their 
dynamics can be modeled by means of effective field theories, allowing a coherent inter- 
pretation of the experimental results. Although very general, this framework can not 
accommodate all phenomena involving elementary particles; in particular, those leading 
to irreversibility and dissipation are clearly excluded. Indeed, a more general treatment is 
needed to properly describe these effects: it can be physically motivated in the framework 
of open quantum systems. [1-3] 

These systems can be thought of as being subsystems in interaction with large en- 
vironments. The time evolution of the total system is unitary and follows the rules of 
ordinary quantum mechanics; nevertheless, the dynamics of the subsystem alone, obtained 
by eliminating the environment degrees of freedom, shows in general irreversibility and 
decoherence. 

When there are no initial correlations between subsystem and environment and their 
mutual interaction is weak, the subdynamics can be described in a mathematically precise 
way in terms of quantum dynamical semigroups. These are linear evolution maps satisfying 
general properties that assure the consistent physical interpretation of the dynamics: they 
include the condition of entropy increase (irreversibility), forward in time composition law 
(semigroup property), complete positivity. This framework is very general and can be 
applied to model irreversibility and dissipation in very different physical situations; [1-11] 
in particular, it can be used to study the evolution of elementary particle systems, treated 
now as open systems. [12-14, 15-19] 

The possibility that decoherence phenomena might affect the physics of elementary 
particles is supported by recent studies on the fundamental dynamics of extended ob- 
jects (strings and branes);[20] indeed, time evolutions described by quantum dynamical 
semigroups can be the result of the interaction with a gas of quanta obeying infinite statis- 
tics (e.g. a gas of D0-branes).[21] In other terms, the dynamics of fundamental extended 
objects could effectively generate at low energies a weakly coupled environment. 

Similar phenomena have also been described in the framework of quantum gravity: due 
to the quantum fluctuation of the gravitational field and the appearance of virtual black 
holes, space-time loses its continuum aspect at distances of the order of Planck's scale and 
assumes a foam like behaviour. [22] As a consequence, new, non-standard phenomena can 
arise, leading to loss of quantum coherence. [23-28] 

Unfortunately, our present knowledge of string theory does not allow to estimate pre- 
cisely the magnitude of the non-standard, dissipative effects induced on elementary particle 
systems; they are nevertheless expected to be very small, being suppressed by at least one 
inverse power of the Planck mass, as rough dimensional analysis suggests. In spite of this, 
the new effects can affect interference phenomena and turn out to be in the reach of future, 
planned experiments. Indeed, detailed investigations of neutral meson systems, neutron in- 
terferometry and photon propagation using quantum dynamical semigroups have already 
allowed deriving order of magnitude limits on some of the phenomenological constants 
parametrizing the new effects, using available experimental data. [29-31, 17, 19] 

In the present work, we shall discuss in detail how non-standard, dissipative phenom- 
ena can affect neutrino propagation, and in particular neutrino oscillations. We shall limit 
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our considerations to the oscillations of two species of neutrinos; in this case, the possible 
dissipative effects can be described in terms of six phenomenological parameters. A pre- 
liminary investigation, limited to vacuum oscillations, has been reported in [18]. There, 
it has been shown that the dissipative phenomena modify the transition probability V 
among the two neutrino flavours, introducing in particular exponential damping factors. 
In a simplified situation, limits on one of the dissipative parameters have subsequently 
been obtained using recent SuperKamiokande data. [32] 

In the following, a much more complete discussion will be presented, with detailed 
analysis of oscillation phenomena in presence of irreversibility, both in vacuum and in 
matter. Dissipation affects both situations; in particular, the resonance condition for 
neutrino propagation in matter turns out to be modified, leading to distinctive observable 
effects. Various approximate expressions for the transition probability V will be given: they 
can be useful in fitting experimental data. A discussion on a possible physical mechanism 
that could give origin to the non-standard effects will also be presented, although much of 
the technical analysis will be relegated to the Appendix. 

As a final remark, let us stress that the presence of non-standard, dissipative phenom- 
ena modify neutrino physics in two important aspects. First of all, they give the neutrinos 
an effective mass, so that oscillations are possible even for massless neutrinos. Further, 
contrary to the standard case, the expression of the transition probability V depends in 
general on the CP-violating phase that is present in the mixing matrix for Majorana 
neutrinos. This allows, at least in principle, to distinguish between Dirac vs Majorana 
neutrinos in oscillation experiments. We find this possibility as one of the most intriguing 
outcome of our investigation. 



2. NEUTRINOS AS OPEN QUANTUM SYSTEMS 

The familiar description of neutrino oscillations involves the study of the evolution of 
neutrinos created in a given flavour by the weak interactions and subsequently detected at 
a later time. The travelling neutrinos are usually assumed to be ultrarelativistic, so that 
the analysis of the transition probability for the original tagged neutrinos to be found in 
a different flavour can be performed using an effective description. [33-37] 

For sake of simplicity, in the following we shall limit our considerations to the mixing 
of two neutrino species. ^ In this case, the neutrino system can be effectively modeled 
by means of a two-dimensional Hilbert space; the two neutrino mass eigenstates will be 
henceforth fixed as basis in this space. In presence of dissipation, the physical neutrino 
states can not be described in terms of elements of the Hilbert space: a more general 
formalism is needed that makes use of density matrices. These are hermitian, positive 
operators (i.e. with non- negative eigenvalues), normalized to have unit trace. 

t The discussion can be generalized to the case of three or more neutrinos; however, the 
explicit formulas for the transition probabilities would become much more involved and 
the discussion less transparent. 
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With respect to the fixed basis, the two flavour states, that we shall conventionally 
call v e and z/ M , are represented by the following 2x2 matrices: 

/ cos 2 # e-^cosfl sin0\ , , 

Pue ~ cos9 sine sin 2 e ) ' ^ Aa) 

( sin 2 # -e-^cosfl sin0\ . , 

pu » ~ y-e 1 * cose sine cos 2 e J- 1 pu ^ [ZAb) 

where e is the "vacuum" mixing angle, while the additional phase ip can be nonvanishing 
for neutrinos of Majorana type. That this extra phase can not be eliminated by a simple 
basis redefinition is a well-known consequence of the reality condition for the Majorana 
neutrinos, and, at least in principle, its presence can be experimentally probed. [38, 39] 
Nevertheless, in the usual approach, this can not happen via the analysis of oscillation 
phenomena alone. [40] As we shall see, the situation is different in presence of dissipative 
effects, so that in the following we shall keep <p nonvanishing unless explicitly stated. 

As explained in the introductory remarks, the evolution in time of any neutrino state 
p will be described by means of linear maps, F t : p(0) i— > p(t), that generalize the stan- 
dard quantum mechanics unitary evolution. Not all generalized maps T t turn out to 
be physically acceptable: they need to satisfy very general physical requirements. First 
of all the maps F t should transform neutrino states into neutrino states, and therefore 
should map any initial density matrix p(0) into a density matrix p(t) = r t [p(0)], for 
any t. Furthermore, they should have the property of obeying the semigroup composi- 
tion law, r t [p(t')] = p(t + £'), for t, t' > 0, of increasing the (von Neumann) entropy, 
S = —Tr[p(t) lnp(t)], of being completely positive. 

It has been proved long ago that evolution maps T t satisfying these properties are 
generated by equations of the following form: [1-3] 

^ = -ztfeff pit) + ip(t) H eS + L[p(t)] . (2.2a) 

The first two pieces in the r.h.s. represent the standard quantum mechanical contributions: 
they give rise to the traditional description of neutrino oscillations in terms of the effective 
(time-independent) hamiltonian -ff e ff- We shall neglect effects due to possible neutrino 
instability: H e g can then be taken to be hermitian. The additional piece L[p] is a linear 
map that encodes possible dissipative, non-standard effects; it can be written as: 

L [p\ = -\Y.{ A ) A iP+p A ) A o) +12 A iP A h ( 2 - 26 ) 

3 3 

where the operators Aj must be such that AjAj is a well-defined 2x2 matrix (entropy 
increase can be easily implemented by taking the Aj to be hermitian). In absence of it, 
pure states (i.e. states of the form \tf;) (tfj\) would be transformed by It into pure states; 
only when the extra piece L[p] is also present, pit) becomes less ordered in time due to 
a mixing-enhancing mechanism: it produces irreversibility and possible loss of quantum 
coherence. 
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In the case of two neutrino flavours, L[p] can be fully parametrized in terms of six, real 
phenomenological constants, a, 6, c, a, /3, and 7, with a, a and 7 non-negative, satisfying 
the following inequalities: [1, 15, 16] 

2R = a + -f - a > , = tfS 

25=a+7-a>0, V = RT 

2T = a + a-7>0, Z = ST 

X = RST - 2 6c/? - #/3 2 - 5c 2 - T6 2 > . 

They are direct consequence of the property of complete positivity. In order for the 2x2 
matrix p(t) to represent a neutrino state, its eigenvalues should be positive for any time 
t; this is crucial for the physical consistency of the whole formalism: the eigenvalues of 
pit) are in fact interpreted as probabilities. The property of complete positivity precisely 
assures that this holds true in any possible condition. (For a complete discussion, see [41].) 

The one-parameter family of finite evolution r t generated by (2.2) are called quantum 
dynamical semigroups; they will be the basis of the phenomenological treatment of the 
dissipative effects in the neutrino system. 

The description of irreversible, non-standard phenomena by means of equations of 
the form (2.2) is actually very general and can be applied to the study of very differ- 
ent physical systems. Originally developed in the framework of quantum optics, [5-7] it 
has also been successfully used in the analysis of statistical models, [1-3] the interaction of 
a microsystem with a measuring apparatus, [8-11] the study of dissipative effects in sys- 
tems involving elementary particles, in particular neutral mesons. [15, 16, 29-31] Although 
essentially phenomenological in nature, all these analysis can be supported by physical 
considerations. 

A general picture in which the quantum dynamical semigroup description of dissipative 
effects naturally emerges is provided by open systems, i.e. by systems in weak interactions 
with a large environment. In the case of elementary particles, these effects are likely to 
originate from the fundamental dynamics of strings or branes, which is in general rather 
complex. Nevertheless, an effective description of the environment that encodes some of 
the properties of the underlying fundamental dynamics turns out to be adequate for a 
more physical discussion of evolution equations of type (2.2). 

Quite in general, the total hamiltonian of a system S in interaction with an environ- 
ment £ can be decomposed as 

H tat = H®l + l®H £ + gH' , (2.4) 

where H is the system hamiltonian in absence of £, while Hs drives the internal dynamics 
of the environment. The interaction between S and £ is described by H' , with g a small, 
dimensionless coupling constant. 

In many instances, the initial state of the total system S + £ can be taken to be in 
factorized form: p tot = p® pe- This is surely justified in the case of the neutrino system: 
since the mechanism of neutrino production is different from the one responsible for the 
dissipative effects, system and environment are surely uncorrelated at the moment of the 



b 2 > , 
- c 2 > , 
p 2 >0, 



(2.3) 
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emission, t Then, the time evolution of the state p of the system S can be obtained by 
tracing over the environment degrees of freedom: 



p = p(0) h- pit) = Tr £ 



-iH t . 



(p®Pe) 



iH tot t 



(2.5) 



In general, the resulting map p(0) h- > p(t) turns out to be rather involved, developing 
non-linearity and memory effects. Nevertheless, when the interaction between S and £ is 
weak, an evolution equation for p(t) local in time naturally emerges. The technical details 
are presented in the Appendix. As discussed there, the environment can be modeled as 
a gas of quanta, obeying infinite statistics; this description is in line with the idea that 
the dissipative effects originate from the low energy string dynamics at a fundamental 
scale Mp (e.g. Planck's mass). Then, in the weak coupling limit, i.e. when the coupling 
constant g becomes very small, the resulting dynamical equation for the subsystem state 
p(t) turns out to be precisely of the form (2. 2). [1-3, 21] 

This result allows a rough estimate of the magnitude of the effects produced by the 
non-standard piece L[p\: they should be proportional to powers of the typical energy of 
the system S, while suppressed by inverse powers of the characteristic energy scale of S. 
In the case of the neutrino system, these effects should be very small, since the typical 
energy scale of the environment can be assimilated to the fundamental scale Mp. For any 
fixed neutrino source and observational conditions, an upper bound on the magnitude of 
the effects induced by L[p] can be evaluated to be of order E 2 /M F: where E is the average 
neutrino energy. 

As a further outcome of the weak coupling limit procedure, the hamiltonian part of 
the evolution equation for p(t) gets modified by the presence of the environment. Indeed, 
the effective hamiltonian H e ^ in (2.2) does not coincide in general with the starting system 
hamiltonian H in (2.4): suitable dissipative contributions to H, generated by the interac- 
tion H' , need to be taken into account. [1-3, 21] As we shall see in the following, this fact 
has interesting consequences in neutrino physics: one can have oscillations among different 
flavours induced by dissipative effects even for massless (or mass-degenerate) neutrinos. 
In other words, originally massless neutrinos can get an effective non-zero mass via the 
interaction with the environment. 



t Even in presence of an initially correlated total system S + £, the factorized approx- 
imation becomes a very good approximation when the short-time correlations have died 
out. [4] 
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3. QUANTUM DYNAMICAL SEMIGROUPS AND NEUTRINO 
OSCILLATIONS 



In the case of the neutrino system, much of the considerations and discussions of the 
previous section about the evolution equation (2.2) can be made more transparent and 
explicit. In particular, both for the effective hamiltonian H e ^ and for the extra piece L[p] 
simple expressions can be given. 

We shall be as general as possible and include in our discussion effects due to the 
propagation of neutrinos in a medium made of ordinary matter. Because of the interactions 
of the neutrinos with the particles in the medium, an effective potential can be generated, 
that has different effects for different flavours. In the case of ordinary matter, the electron 
neutrinos interact with the electrons in the medium, so that their average energy effectively 
receive an extra contribution A = \/2Gf n e with respect to the energy of the muon 
neutrinos (Gp is the Fermi constant, while n e represents the electron number density in 
the medium). [42, 33-37] In the ordinary case, this contribution can significally change the 
oscillation pattern between v e and states (the so-called MSW effect). [43, 44] As we 
shall see, this phenomenon can be substantially modified by the presence of non-standard, 
dissipative effects. 

In the basis introduced in the previous section, the 2x2 matrix representing the 
effective hamiltonian can be taken to be of the form: 



In the first piece, E represent the average neutrino energy, while ujq = Am 2 /4E takes into 
account the square mass difference Am 2 of the two mass eigenstates; these are the usual 
contributions that give rise to the standard oscillation pattern in vacuum. The extra real 
parameters ui, u)2 and us are the consequence of the interaction with the environment; as 
explained in the previous section (and discussed in detail in the Appendix), they represent 
the contribution of the dissipative phenomena to the system hamiltonian. 

Both ujq and u\, uj 2 , cu 3 contribute to the level splitting u> = [(uj + uj 3 ) 2 + uj 2 +W 2 ] 1 / 2 
between the two mass eigenstates, so that they all contribute to the oscillation phenomena 
in vacuum. Therefore, even for initially degenerate mass eigenstates, Am 2 = 0, vacuum 
oscillations can occur between the two flavours due to the dissipative effects induced by 
the fundamental dynamics at the large scale M F . Although in general all three parameters 
u>i, UJ2 and UJ3 are non- vanishing, in the following, in order to simplify the treatment, we 
shall assume u\ = 0J2 = 0; this working assumption allows for more manageable formulas, 
while keeping unaffected their physical meaning and implications. ^ 

The final contribution to H e ^ in (3.1) takes into account the interaction of the prop- 
agating neutrinos with ordinary matter; it would be diagonal in the flavour basis (only 
electron neutrinos are affected), but assumes a more complicated expression involving the 
mixing angle 9 and the phase f in the chosen basis. Since the coefficient A is proportional 
to the density of electrons in the mean, for propagation in non-homogeneous matter i7 e ff 

t When Am 2 = 0, this is no longer an assumption: in this case, one can always choose 
to work in a basis for which uii and U2 vanish. 



I — (jjg — CU3 



Ui — IUO2 

E + ujq + uj 3 




A/l + cos2# e-'v sin 29 
2 V e** sin 29 1 - cos 29 



(3.1) 
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will in general be a function of the position of the neutrinos. Nevertheless, one can always 
approximate a non-homogeneous medium by a collection of media, each with a constant 
density, while having different thickness; in view of this, in the following we shall assume 
the parameter A to be a constant (see also the discussion in Sect. 4). 

As mentioned before, although the effective hamiltonian H e ff gets also dissipative 
contributions, only when the additional piece L[p] in the evolution equation (2.2) is non- 
vanishing, irreversibility and mixing enhancing effects are possible. In the present case, 
its explicit expression in terms of the six phenomenological constants a, 6, c, a, /?, and 7 
in (2.3) can be most simply given by expanding the 2x2 matrix p in terms of the Pauli 
matrices Uj, i = 1,2, 3, and the identity uq: 



1 3 

( 3 - 2 ) 



2 

,u=0 



In this way, the linear map L acting on p can be represented by the following, symmetric 
4x4 matrix [L^,] , acting on the 4- vector of components (po, p\, pi, pz)'- 



[V] = " 2 











°\ 





a 


b 


c 





b 


a 


p 


Vo 


c 




1) 



(3.3) 



The form of the evolution equation (2.2) can be further simplified by recalling that it 
is trace preserving. From the initial normalization condition Tr[p(0)] = 1, one immediately 
obtains that the component of p(t) along the identity is equal to one for all times. Then, 
the evolution equation for the remaining three components of p(t) can be rewritten in a 
Scrodinger-like form: 

^\p(t)) = -2H \p(t)) , (3.4) 
where the 3- vector \p) has components (pi, P2, P3), while 

(a b + [i c — v sin ip 

b — p a (3 + v cos ip ) , (3.5) 

c + v sin (p f3 — v cos p> 7 

with 

A A 
p = — cos 29 — u , v = — sin 29 . (3-6) 

The solution of (3.4) involves the formal exponentiation of the matrix H: 

\p(t))=M(t) |p(0)> , M(t)=e- 2nt . (3.7) 

As discussed in [18], expressions for the entries of Ai(t) can always be obtained by 
solving the eigenvalue problem for the 3x3 matrix in (3.5): 

H\v {k) ) = A (fc) |f (fc) > , k = 1,2,3. (3.8) 
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The three eigenvalues X^ 2 \ X^ satisfy the cubic equation: 

X 3 + rX 2 + sX + w=0, 

with real coefficients: 



(3.9) 



r = -(A (1) + A (2) + A (3) ) = -(a + a + 7) , 



(3.10a) 



s ee A^A^ + A^A^ + X^X^ = aa + a 1 + a 1 -b 2 -c 2 -(3 2 + p 2 + u 2 , (3.106) 



w ee -A (1) A (2) A (3) = a(/3 2 - z/ 2 cos 2 + a(c 2 - v 2 sin 2 + 7 (o 2 - ^ 2 ) 
— aa7 — 2 6c/3 — 6z/ 2 sin 2(p — 2p,v(j3 sin ip + c cos </?) 



(3.10c) 



The solutions are either real, or one is real and the remaining two are complex conjugate, 
according to the sign of the associated discriminant: D = p 3 + g 2 , p = s/3 — (^/3) 2 , 
q = (r/3) 3 — rs/6 + w/2 (degenerate, real solutions occur when D = 0).[45] Then, recalling 
that the matrix H itself satisfy the equation (3.9), one can derive the following expression 
for the entries of A4(t): 



-2A (fc) t 



fc=l 



]2 +rA (fc) +s ) ( J y+ ( A (fc) +r ) W . j . +W I 



13 



3[A( fc )] 2 + 2rA« + s 



z,j = l,2,3. 
(3.11) 

Although rather formal, this formula allows a general discussion on the behaviour 
of M.(t). For fj, = v = 0, thanks to the inequalities in (2.3), the matrix H results real, 
symmetric and non-negative: its eigenvalues are all real and non-negative. Only when 
\fi\ and \u\ are sufficiently large, complex eigenvalues may appear, although with a non- 
negative real part, since in general the evolution generated by (2.2) is bounded for any 
£.[46] In this case an oscillatory behaviour is possible, while for small n, u, the damping 
terms prevail and dissipation is the dominant phenomena. 

In particular, since generically detTC = —w 7^ 0, in presence of dissipation the real part 
of A^ 1 ), A^ 2 ), A^ 3 ) are all strictly positive; therefore M(t) asymptotically vanish for large 
enough times. t This has clearly dramatic consequences in the study of neutrino flavour 
transitions. 

Let us assume that at t = the neutrinos are generated to be of type v e . In the 
formalism of density matrices, the probability of having a transition into neutrinos of type 
at time t is given by: 



1 + EpiX m *&) 



i,3=l 



(3.12) 



t In presence of vanishing eigenvalues, this decoherence effect is only partial; [18] how- 
ever, note that having detH = requires a unnatural fine-tuning among the parameters 
in (3.10c). 
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where p Va if) is the solution of (2.2) with the initial condition given by the matrix p Ue (0) = 
p„ e , while p^ e , pi , i,j = 1, 2, 3 are the components of the 3- vectors \p Ve ), \pu^) correspond- 
ing to the density matrices in (2.1). Using the explicit expressions for these components, 
one finds: 

r, e ^(t) =^|l- cos 2 29 M 33 (t) 

-sin 2 26> Mu(t) cos 2 + M 2 2(t) sin 2 <p + (m 12(f) + M 2 i(f) s j simpcostp 
cos29 sin26 \(Mi 3 (t) + «M 3 i(*)) cos <p + (M 2 3(t) + M 3 2(tf) sinp 



(3.13) 

One of the interesting features of this formula is its explicit dependence on the phase <p; 
in presence of dissipative effects, it is therefore possible, at least in principle, to distinguish 
between Dirac and Majorana neutrinos by studying the oscillation pattern in (3.13). This 
peculiarity disappears when the non-standard, dissipative pieces in (2.2) are absent; indeed, 
in that case, one has: 

Mij(t) = Stj Htj H 2 ^ *j ^7 J = 1 5 2, 3 , (3.14) 



where 7i is now as in (3.5) with a, 6, c, a, and 7 all equal to zero, while % = y p 2 + v 2 , 
and (3.13) reduces to the well-known standard expression for the oscillation probability in 
an homogeneous medium: [42, 33-37] 

Vl%^ (t) = sin 2 29 M sin 2 iv M t , sin 2 26> M = • (3.15) 



(^-cos2#) + sin 2 20 



Another distinctive characteristic of the transition probability in presence of dissi- 
pation given in (3.13) is its asymptotic behaviour for large times, which turns out to be 
independent from the mixing angle 9, the phase ip and the matter coefficient A: 

V„ e ^(t) ~ \. (3.16) 

This result is a direct consequence of the vanishing of the matrix M.(f) in (3.7). Never- 
theless, as discussed below, the regime of validity of this asymptotic limit can seldom be 
reached in practical experimental conditions. 
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4. TRANSITION PROBABILITY IN MATTER 

The general expression (3.13) for the transition probability is rather involved and it is 
not particularly useful for studying in more detail its physical properties. Therefore, in the 
present and following sections we shall discuss various approximations in which V Vf ,^ v (t) 
assumes a more manageable form. These simplified expressions, besides being appropriate 
for theoretical analysis, could also be used to fit actual experimental data. [32, 47] 

As already mentioned in Sect. 2, the values of the constants a, b, c, a, j3, and 7 
parametrizing the non-standard effects, are expected to be very small, with an upper bound 
of order E 2 /M F ~ 1(T 19 GeV, for E ~ 1 GeV and for M F the Planck mass. Nevertheless, 
this estimate is not far from the values that the standard oscillation parameter ujq = 
Am 2 /4E assumes for typical neutrino sources. Indeed, the ratio of a, b, c, a, (3, and 7 with 
ujo can be evaluated to be at most of order 10~ 10 E 3 / Am 2 , with E expressed in MeV and the 
neutrino mass difference Am 2 in eV 2 ; this ratio turns out to be about 10 2 for atmospheric 
neutrinos, of order one for solar neutrinos, while for accelerator neutrinos it can be as small 
as 10 -2 . Therefore, the effects induced by dissipation can interfere with those producing 
oscillations via a non- vanishing ujq, resulting in observable modifications of the oscillation 
pattern. Present and, most likely, future dedicated neutrino experiments should be able 
to detect these modifications, or at least put stringent limits on the magnitude of the 
non-standard phenomena. 

Let us first consider the case in which the dissipative parameters a, b, c, a, j3, and 7 are 
of the same order or larger than the remaining constants in (3.1). In this case a very useful 
approximation is to assume a = a = 7 and c = 0, conditions perfectly compatible with 
the inequalities (2.3), provided a 2 > b 2 + f3 2 . For simplicity, we further assume the extra 
phase (p to be vanishingly small. A manageable expression for the transition probability 
can then be derived: 

e~ 2at sin 2 (0 M t) , (4.1) 
where 

n M =[^ 2 + u 2 -b 2 -(3 2 } 1/ \ (4.2) 
9 = u;sm20, (3 = (3 cos26 + b sin2# . (4.3) 

The oscillating behaviour in (4.1) depends on the magnitude of the combination \i 2 + is 2 = 
(A/2— 00 cos2#) 2 +uj 2 sin 2 20 with respect to b 2 +(3 2 ; in regions for which b 2 + (3 2 > [i 2j rv 2 , 
the frequency Qm becomes purely imaginary and V Va ^ v (t) contains only exponential 
terms. Anyway, the a-depending damping terms in (4.1) dominate for large times, and 
the asymptotic limit (3.16) is thus recovered. 

In absence of dissipation, a = b = f3 = 0, the factor in front of the sine term in 
(4.1) becomes parametrized as in (3.15), with a modified mixing angle 9m- When the 
matter parameter A is close to Ar = 2uj cos 20, the transition probability gets enhanced, 
and oscillations between the two neutrino species is possible even when the original mixing 
angle is small. This phenomenon is at the root of the so-called MSW effect. [43, 44, 42] 

In presence of dissipation however, the physical consequences of this effect are in 
general much more modest. In this case, one can not parametrize V v ^ v (t) in terms of a 
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modified mixing angle; in spite of this, the expression in (4.1) as a function of A, at fixed 
time, has a critical point for A = Ar. This point is a maximum for v 2 > /3 2 , and indeed as 
A approaches Ar an enhancement in V Ve ^v (t) occurs: (z/ 2 — /3 2 )/0 2 kf > 1; however, the 
exponentially damping factors in (4.1) greatly reduce in practice its effectiveness. Further, 
when v 2 < /3 2 , the probability Vv e ->v^(t) in (4.1) is maximally suppressed at the critical 
point: it is dominated by the damping factors. 

This discussion might appear spoiled by the initial assumption of a constant matter 
parameter A: the occurrence of the MSW effect requires a medium with a (slowly) varying 
density. As already pointed out, the assumption of a constant A is not really a limitation: 
one can always approximate, with arbitrary accuracy, the travelling of neutrinos through 
varying density matter as the propagation in a series of media with different constant 
densities and different thickness. The total time-evolution will be given by the composition 
of the evolutions in the various matter slices, so that the matrix M. in (3.7) becomes: 

M(t) =M n (t n )-- ■ M 2 (t 2 ) Miih) , t = t 1 +t 2 + ... + t n , (4.4) 

where, ti, t 2 , . . . ,t n are the total times spent by the neutrinos in the various media, while 
Mi, i = 1, . . . , n are the corresponding propagation matrices. 

As an example, let us consider the case of an initial electron neutrino travelling for 
a time t\ into a medium with matter parameter A, that is then detected in vacuum at a 
later time t = t\ + 1 2 ; this situation can roughly represent a solar neutrino model. Using 
(4.4), the probability of detecting the original v e as a muon neutrino is given by: 



_ e -a(ti+t 2 ) 



j> 2 -/3 2 

n 2 



4fi 2 - AM + 26 + 2a; cos 20) . 2// _ . . 2/r ^ . 

M 



o o 2 

-— - sin(20 M ti) sin(2O ^) - 2sin 2 (O t2) - 2-^- sin 2 (fi M ^ 
il M \l M 

(4-5) 

where Oo = \f^ 2 — b 2 — /3 2 , while Om is as in (4.2). One can check that, with the appro- 
priate choice of parameters (see the previous discussion), the probability V Ve ^ v ^{t) can 
indeed get an enhancement for A close to the critical point; the effect is however modest 
and further suppressed by the damping factors. Nevertheless, with the appropriate values 
of A, t\ and t 2 , the expression (4.5) can be used to fit solar neutrino data; the total flight 
time t is large, so that a good sensitivity at least on the dissipative parameter a is surely 
attainable. 

A different approximation of the full expression (3.13) for the transition probability 
Vv e ^v can be obtained when the dissipative parameters a, 6, c, a, /3, and 7, can be 
considered small with respect to the level-splitting term uj; as mentioned before, this 
typically occur for neutrino beams generated at accelerators. In this case the additional 
term L[p] in the evolution equation (2.2) can be treated as a perturbation. To first order 
in the small parameters, explicit expressions for the entries of the evolution matrix M. (t) 
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in (3.7) can be easily obtained; then, using (3.12), one finds: 



V 2 \ e -2X lt + e -X 2 t 



uj 



hi 



V 2 



UJ 



M 



COs(2 UJMt) + 



N 



UJ 



M 



where ujm = a/ A* 2 + 1/2 as m the previous section, while 



sin(2 UJMt) 



, (4.6) 



p, = fi cos 29 + v sin 29 = — — uj cos 29 , 
v = v cos 29 — /j sin 29 = uj sin 29 ; 



(4.7a) 
(4.76) 



the parameters Ai, A2 and N contains the dependence on the dissipative constants: 

Ai = [av 2 + 2c\±v + ^h 2 )/uj 2 m , 
A 2 = a + [an 2 - 2c[iv + ^u 2 )uj 2 m , 



(4.8a) 
(4.86) 



v 2 

N = — 
2 



a — a — v(3av + 2cfi)/u' M + 2>vv{a — 7) cos 29 



+ c 



(2fip, 2 v - uj m sin 49 ) /uo 2 M . (4.8c) 



In the expression (4.6), we have reconstructed the exponential factors by consistently 
putting together the terms linear in t. Notice that the result in (4.8) is in agreement 
with the discussion in Sect. 3 concerning the eigenvalues of the matrix H. In this case the 
algebraic equation (3.9) has one real, A^, and two complex conjugate solutions, A^ 2 ' 3 ^ = 
Xr ± iAj; within our approximation, A^ 1 ) = Ai, 2Ai? = A2, so that the first condition in 
(3.10) is satisfied, while the remaining two fix the imaginary part A/. 

The expression (4.6) for the transition probability can be used to fit experimental 
data. With respect to the standard case, it contains three additional parameters, Ai, A 2 
and N, that signal the presence of non-standard phenomena, through the constants a, 
6, c, a, /?, and 7. If at least one of these three parameters is found to be non zero, it 
would clearly signal the presence of dissipative effects in neutrino physics; this is surely 
the most simple experimental check on the generalized evolution equation (2.2) that can 
be performed with accelerator neutrino beams. 



5. ADIABATIC APPROXIMATION 

For time evolutions with a semigroup composition law, the appropriate way to follow 
neutrino propagation in a medium is through the successive applications of the finite 
evolution matrices Ai(t) to the initial state |p(0)), as shown in (4.4). In more traditional 
approaches, one usually adopts a different approximation, based on the assumption of 
(adiabatic) slowly varying matter density. This approximation can be easily discussed also 
in the framework of density matrices and quantum dynamical semigroups. 
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Let us consider the case of a neutrino, created at t = in matter of high density (i.e. 
in the core of the sun), propagating towards regions of smaller density. In this case, the 
effective hamiltonian (3.1) is no longer constant, and the propagating matrix -M(£) involves 
a time-ordered exponentiation of 7i. Nevertheless, at any instant of time, the 3x3 matrix 
Tt can be diagonalized by a similarity transformation: 

H = TVT~ X . (5.1) 



Using this decomposition in (3.4), one can derive the evolution equation for the transformed 
3- vector, \p) = T |p); explicitly, one finds 



>«)> = "2 



at 



\m) ■ (5-2) 



The adiabatic approximation amounts to neglecting the last term in this equation; this is 
justified when the matter density parameter A is slowly varying. In this approximation, 
the neutrino state essentially evolves in time as an eigenstate of 7i. 

In order to make the discussion more explicit, as in the previous section we shall take 
a = a = 7 and c = 0, while neglecting the extra phase <p and the dissipative contributions 
ui, u)2 to H. With these choices, one has: 

a 

V=-2\ a + m M | , (5.3) 

a — iQm 

with Qm as i n (4-2), while the transformation matrix T takes the form: 

1 I V^((3 + u) b + /j b + /j 

ifl M —i£Im I • (5.4) 



The entries of T can be parametrized in terms of two real variables £ and £ and an angle 
4>, that could be complex: 



M 



b±p 



M 



■■ ± (e ±c cos 29 sin 20 + e ± ^ sin 29 cos 2<^j , (5.5a) 
±(^e TC sin2^sin2(/)-e T? cos2^cos2(/)j ; (5.56) 



for later convenience, an explicit dependence on the mixing angle 9 has been extracted 
from the entries of T. 

Using a compact vector notation, the transition probability V v< ,^, v ^ in (3.12) can now 
be written as 

V„^{t) = \{l + ( P ^\T r M- D {tf,t i )-T- 1 \p Vu )} , (5.6) 



14 



where 

Mv{t f ,U)=e- 2at [ e -*£ t ) , (5.7) 

2i (ItQm(t) 



e 



while Ti and Tf are the matrices that diagonalize H. at the initial time ti = and final 
time tf = t; they can be written as in (5.4), with parameters Q, 0, and £/, Q, cf>f, 
respectively. Explicitly, one finds 



V„ e ^(t) = ^|l-e" 2Q!t [e^-^ cos 20/ cos 20 



+ e -(C/-Ci) 

sin20/ sin20j cos ^2 / rfrO M (T)j 



(5.8) 



In absence of dissipation, £j = Q = £/ = Q = 0, a = b = (3 = 0, one recovers the familiar 
expression for the adiabatic transition probability. 

When the adiabatic approximation ceases to be valid, the previous treatment needs 
to be generalized. Indeed, in this case, the neutrino state no longer remains in a spe- 
cific eigenstate of Tt for the whole time evolution; rather, it can mix with the remaining 
eigenstates. In order to take into account this possibility, the expression for the transition 
probability (5.6) needs to be modified: 

Vu^u^t) = i{l + (pvj T f ■ Mv(t f ,t c ) ■ A ■ Mv(t c ,U) ■ T- 1 |p„ e )} , (5.9) 

where t c is the time at which the neutrino crosses the critical region, while A is the 
mixing matrix that encodes the possible hopping between the instantaneous eigenstates 
of the effective hamiltonian. For simplicity, we are not taking into account the possibility 
of hoppings induced by dissipative effects: they can be considered to be negligible with 
respect to the matter induced ones. As a consequence, A can be taken to be the most 
general 3x3 unitary matrix, that preserves appropriate consistent conditions: they assure 
the reality of the transition probability. Taking into account these conditions, A can be 
parametrized in terms of two complex numbers u and v, such that \u\ 2 + \v\ 2 = 1: 

(\u\ 2 — \v | 2 \pluv \pluv\ 
-V2uv u 2 -v 2 . (5.10) 
-\[2uv -v 2 u 2 ) 

The explicit expression for the probability in (5.9) is now rather involved; however, it 
simplifies when neglecting the fast oscillating terms: 

Vu e ^{t)) = ^{l - e~ 2at e € '-fc (1-2 M 2 ) cos 20/ COS20,} . (5.11) 

In practical applications, the interesting case occurs when the neutrinos are generated 
in a medium with very large matter density, cos20i ~ — 1, while detected at a later 
time t in vacuum, £/ ~ 0, 0/ = 9. In this case, one finds: 

V^ u ,(t)) = i(l-e" 2at ) +e- 2at [cos 2 #- M 2 cos2#] . (5.12) 
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This is the most simple form that the transition probability formula takes in presence of 
dissipative and matter effects: with respect to the familiar expression, it contains expo- 
nential damping factors. Taking into account that neutrinos are relativistic, the flight time 
between emission and detection is with very good approximation the same as the distance 
£ between source and detector. One can then use (5.12) to derive a rough order of mag- 
nitude limits on the non-standard parameter a. The best bounds are expected from solar 
neutrinos, where l/£ can be as low as 10~ 27 GeV. 



6. TRANSITION PROBABILITY IN VACUUM 

One of the most interesting properties of the quantum dynamical semigroup approach 
to neutrino propagation is the possibility of probing the nature of the neutrinos by studying 
their oscillation pattern. Indeed, the transition probability V Va -^ v in (3.13) explicitly 
depends on the extra phase <p, which can be non-vanishing for Majorana neutrinos. In this 
section we shall discuss to what extent the phase ip can be extracted from V Ve ^> v . With in 
mind possible applications to atmospheric neutrinos, we shall limit our considerations to 
oscillations in vacuum, and present various explicit formulas for the transition probability 
in different approximations. 

As a working assumption, let us first take the dissipative parameters c and (3 in (3.5) 
to be much smaller than the remaining constants. To lowest order, a closed form for the 
entries of the evolution matrix Ai(t) in (3.8) can then be obtained. Using the general 
formula (3.13), one finds: 

V v ^ (t) = -{l -e- 27t cos 2 26'-e- (a+a)t sin 2 26'[cos(2 0t) 

\B\ 11 (0) 

+ ^ sin(2 0t) cos(0 B + 2v?)J j , 



where B = a — a + 2ib = |£?|e l ^ B and O = a/cu 2 — |-B| 2 /4. In this case the dependence 
on the phase <p> is very mild, and can not be extracted by studying V v ^ v alone: an 
independent determination of the combination B is necessary. 

The situation is even worse when 7 = 0; in this case, the inequalities (2.3) automat- 
ically guarantee c = f3 = and further impose b = 0, a = a. In this case, ip completely 
disappears from the expression of the transition probability, that reduces to: [18] 

Vu^up (t) = X - sin 2 29 [l - e~ 2at cos(2 ut)j . (6.2) 

In view of this, analysis of the experimental data based on (6.2) along the lines of Ref.[32] 
is totally insensitive to cp; a fit with more than one non-vanishing dissipative parameters 
is in general needed, although this condition is certainly not enough, as shown by (6.1). 
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In this respect, a more interesting situation occurs when c = and a = a = 7, as 
considered in the previous sections. All the entries of the evolution matrix Ai(t) are now 
non-vanishing, and the transition probability can be written as: 



V v .^{t) = \{ l - z~ 2at ( 1 + ^T sin 2 (fi t)) cos 2 29 



( 2/3 2 b 

+ sin 2 26* I cos(2O ^) — 7^2 sm2 (^o^) cos 2 (p — — sin(20o^) sin2(/? 
V "0 

+ — cos 49 sin(Oo^) ( tt - sin(O t) cos <p — cos(O t) sin <p 

"0 \"0 

(6-3) 

where fi = \Joo 2 — b 2 — j3 2 as before. In this case, the dependence of V Ve ^ v (f) on time 
is significally altered by the presence of a no n- vanishing <p. 

Atmospheric neutrino data are the most suitable for an experimental study of (6.3), 
since in this case the time dependence can actually be probed. Nevertheless, it should 
be stressed that the expression in (6.3) contains four additional parameters besides the 
standard ones, u and 9, so that the fitting procedure might turn out to be difficult in 
practice. In order to simplify the analysis, one can further assume one of the two non- 
standard parameters b or (3 to be zero, but not both: here again, when a is the only non- 
vanishing dissipative parameter, the dependence on <p in V Ue ->v (t) disappears. Despite 
these difficulties, the amount of data on atmospheric neutrinos is constantly growing, so 
that at least some information on the presence of a non-vanishing <p, together with some 
of the dissipative parameters, will surely be attainable in the near future. 



7. DISCUSSION 

The study of open systems by means of quantum dynamical semigroups offers a physi- 
cally consistent, general approach to the discussion of phenomena leading to irreversibility 
and dissipation. When this formalism is applied to the analysis of the propagation of 
neutrinos, both in vacuum and in matter, it gives precise predictions on the pattern of 
oscillation phenomena: the new, non-standard effects manifest themselves through a set 
of phenomenological parameters, ui, U2, uj^ (the hamiltonian ones), a, 6, c, a, [3, and 
7 (the purely dissipative ones). Their presence allows oscillating phenomena even for 
mass-degenerate neutrinos, accompanied by possible CP-violating effects. Further, these 
predictions can be experimentally probed. Indeed, fits of experimental data along the 
lines discussed in [32] can be repeated for the more complete expressions of the transition 
probability V Va ^ v (t) presented in the previous sections. 

For sure, the fitting procedure will be more difficult and uncertain than in the standard 
case, due to the presence of more unknown parameters. Furthermore, one has to take into 
account that the effects induced by the presence of the new parameters are expected 
to be small. By viewing the neutrinos as an open system in weak interaction with an 
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environment generated by a fundamental "stringy" dynamics, the effects of irreversibility 
and dissipation can be roughly estimated to be proportional to the square of the average 
neutrino energy, divided by the characteristic energy scale of the environment. Assimilating 
this scale to the Planck mass produces estimates of order 10 -27 GeV for solar neutrinos, 
while larger values are expected for more energetic neutrinos. Despite these difficulties, 
a new generation of dedicated neutrino experiments are presently collecting data or will 
shortly start construction, so that stringent bounds on the dissipative effects can surely be 
expected in the future. 

The above estimate on the magnitude of the non-standard, dissipative effects is based 
on very general and physically motivated considerations about open systems; therefore, it 
is rather robust and quite independent from the details of the microscopic, fundamental 
dynamics responsible for the interaction between the neutrino subsystem and the envi- 
ronment. Nevertheless, it has been questioned on the basis of a formal similarity of a 
particular, simplified version of the evolution equation (2.2) with those describing the 
phenomenon of the so-called dynamical reduction of the wave-packet. [48] The analogy is 
rather superficial: the physical process leading to dissipation and the reduction process are 
quite distinct and act at different energy scales. Furthermore, as a more complete analysis 
would reveal, quantum dynamical semigroups generated by equations of the form (2.2) are 
unable to properly describe dynamical reduction processes. [9] 

A different criticism on the use of quantum dynamical semigroups for the description 
of dissipative effects advocates the use of non-linear evolution equations. [49] Once more 
the general theory of open systems offers a clarifying discussion on this point (for further 
technical details, see the Appendix). 

As pointed out in Sect. 2, the dynamics of a small system S in interaction with a 
large environment £ is in general very complex and can not be described by means of 
evolution equations that are linear in time: possible initial correlations and the continuos 
exchange of energy as well as entropy between S and £ produces memory effects and 
non-linear phenomena. Nevertheless, when the typical time scale in the evolution of the 
subsystem S is much larger than the characteristic time correlations in the environment, 
the subdynamics simplifies and a mathematically precise description in terms of quantum 
dynamical semigroups naturally emerges. [1-3, 21] 

This limiting procedure is general and can be applied to all physical situations for 
which the interaction between S and £ is weak and for not too-short times, so that the 
non-linear disturbances due to possible initial correlations have died out. [4] These are pre- 
cisely the conditions that are expected to be fulfilled in the case of neutrino systems: the 
characteristic time correlations in the environment, induced by the fundamental (gravita- 
tional or stringy) dynamics, is certainly much smaller than the neutrino propagation time, 
while the interaction between neutrinos and environment is for sure weak. 

The description of neutrino propagation in terms of quantum dynamical semigroups 
automatically guarantees the fulfillment of basic physical requirements, as forward in time 
composition, entropy increase (irreversibility) and complete positivity. This is a clear 
advantage over alternative formulations. Based on ideas originally presented in [23], gen- 
eralized dynamics for the neutrino system incorporating some of these properties have 
been discussed before (see [50-52]). However, those dynamics do not satisfy the condition 
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of complete positivity; as already mentioned, this could lead to serious inconsistencies that 
can be avoided in all situations only by adopting evolutions of the form (2.2). [41] 

Non-linear dynamics in the description of the neutrino system naturally emerge when 
the requirement of weak coupling between neutrinos and environment is not satisfied. This 
typically happens in extreme conditions, as those found in the core of a supernova or the 
early universe; more conventional dissipative phenomena then arise due to the scattering 
and absorption processes in the medium. [53-55] In order to properly deal with these situa- 
tions, a second-quantized, field-theoretical formalism has been constructed, using specific 
effective interaction hamiltonians as starting point. Although derived using techniques 
similar to the ones described before, the resulting kinetic evolution equations are quite 
distinct from (2.2); they give rise to decoherence effects that modifies the pattern of neu- 
trino oscillations in a very different way with respect to the expressions discussed in the 
previous sections. Nevertheless, also in this case the condition of complete positivity needs 
to be satisfied for consistency, and this requirement might produce further constraints on 
the modified dynamics. 

Decoherence effects in neutrino physics have been further discussed in connection 
with the uncertainties in the emission and detection processes. [56] By smearing the famil- 
iar expression for the transition probability over energy and time (or position) with an 
appropriate Gaussian distribution, an exponential damping factor is generated, so that 
the resulting expression for the averaged probability looks similar to the one presented in 
(6.2). The analogy is once more only superficial, since the transition probability (6.2) (and 
more generally the expression in (3.13)) has an explicit time (position) dependence that 
can not be reproduced via a Gaussian average. Further, the physical mechanisms leading 
to the modified probability expressions are clearly different, the detector "noise" in one 
case, a fundamental dynamics in the other; in turn, this leads to a different dependence of 
the damping factors on the average neutrino energy. 

As mentioned before, quantum dynamical semigroups can be employed to model a 
large variety of physical situations. It is not a surprise that they have been used to 
study the effects of density waves in the propagation of neutrinos in fluctuating media, 
in particular, in the interior of the sun; these phenomena are also described by equations 
of the form (2.2) and induce modifications on the neutrino oscillation probabilities. [57] 
However, it should be stressed that these density fluctuations have their origin in the 
dynamics of the sun and operate at energy scales quite different from the Planck mass; 
therefore, they can be easily isolated from the dissipative effects discussed before, that are 
not expected to be influenced by long-range phenomena. 

As a final remark, let us point out that several unconventional phenomena affecting 
neutrino propagation have been discussed in the literature; they include: neutrino decay, 
flavour changing neutral currents, violation of Lorentz invariance or the equivalent principle 
(e.g. see [58-64]). All these phenomena lead to modifications of the standard oscillation 
pattern; however, the resulting transition probability V Ue ^ u ^(t) has a dependence on time 
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(or pathlength) and neutrino energy that differ from the one discussed in the previous 
sections. 

Indeed, the dependence of the observable V Vt ,^ v it) on the phenomenological param- 
eters a, 6, c, a, (3 and 7 is very distinctive of the presence of dissipative phenomena and 
can not be mimicked by other unconventional mechanisms. This further strengthens the 
possibility of identifying the dissipative contributions from the analysis of experimental 
data, quite independently from other effects. 



APPENDIX: THE WEAK COUPLING LIMIT 

For a neutrino system in interaction with an environment S, the total Hamiltonian 
can be decomposed as in (2.4): 

fltot = H®l + l®H £ + gH' , (Al) 

where, neglecting for simplicity matter effects, the system hamiltonian H can be taken to 
be of the form 

Assuming no initial correlations between the two systems, the evolution in time of the 
neutrino state p(t) follows the general rule (2.5): 

p p{t) = Tr £ \ e - iHtott (p <g> ps) e iHtott ] , (A3) 

This evolution map is in general very complicated, developing irreversibility and memory 
effects. However, it simplifies when the interaction between the neutrino subsystem and 
the environment is weak. 

There are essentially two different ways of implementing in practice this condition: [1-3] 
they correspond to the two ways of making the ratio t/t£ large; here, r is the typical 
variation time of pit), while T£ represents the typical decay time of the correlations in the 
environment. Only when r ^> Tg, one expects the memory effects in (A2) to be negligible, 
and a local in time evolution for the state p(t) to be valid. 

When ts becomes small, while r remains finite, one speaks of "singular coupling limit" , 
since the typical time-correlations of the environment approach a 5-function. In the other 
case, it is r that become large, while T£ remains finite: one then works in the framework 
of the "weak coupling limit". In practice, this is obtained by suitable rescaling the time 
variable, t — > t/g 2 , and by sending the coupling constant g to zero (van Hove limit). [1-3] 

The choice between these two limiting procedures clearly depends on specific physical 
considerations about the system under study. In the case of the neutrino system, both 
limits appear physically acceptable. Since the singular coupling limit has been presented 
elsewhere, [21] in the following we shall concentrate on the discussion of the weak coupling 
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limit. In this case, following the steps presented in [21], from the finite time evolution 
(A. 3) one can derive a differential equation for p(t) local in time; it is of the form (2.2a), 



^ = -iH p(t) + ip(t) H + C[p(t)} , (A4) 



with 



C[p] = - lim / ds [ dt Tr g \e lH ^ s \e iHot H' e~ iHot , \H' , p ® p £ ]\e- %H ^ s \ , 
T^oo 2T J_ T J [ L J J 

(A5) 

where i^o represents the limit of H tot for a vanishing coupling constant g. 

The general form of the additional term C[p] in (A. 5) does not actually depend very 
much on the details of the environment dynamics; an effective description that takes into 
account its most fundamental characteristic properties is enough to allow an explicit eval- 
uation of the integrals in (A5). Following the idea that the dissipative effects are low 
energy phenomena that originate from the fundamental gravitational or stringy dynamics 
at some large scale Mp, we shall model the environment as a gas of quanta, obeying infinite 
statistics, in thermodynamic equilibrium at inverse temperature (3p = 

Further, taking into account that the interaction between the neutrino system and the 
environment is weak, we shall assume the interaction hamiltonian H' to be linear both in 
the neutrino and the environment dynamical variables: 

3 

an explicit expression for the environmental operators will be discussed below. 

In order to proceed further, it is convenient to introduce a spectral decomposition, 
and use the auxiliary matrices crjf\ A = — 1, 0, 1,[1] 

4°) = P 1 a fl P 1 +P 2 a^P 2 , a^=P 1 a^P 2 , a^ = P 2 a^P 1 , (A.7) 

with 

*-(;s). *-(!!!)■ <-) 

Then, the limit in (A.5) can be explicitly performed and the result expressed in terms of 
the following 4x4 hermitian matrices: 



a$ =g 2 dt e" 2 ^ (B,(t) B v ) , (A.9a) 

J — oo 

b$=ig 2 {f~dte- 2iXuot (B„(t)B„)- dt e 2 ^ (B^B u (t))} , (AM) 



t For a motivation of this choice in terms of the dynamics of extended objects (D0- 
branes), see [65-67, 21]. 
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involving thermal correlations of the environment operators, 



(B,(t) B v ) = Tr £ [B.it) B v p £ ] , p £ = ^z^) ■ 
Explicitly, one finds: 



C[p] = i[p , H] + L[p] , 



where 



iw = 5 E E< A 



(A) (A) (A) (A) (A) (A)' 

2 a) 'pal - a\ 'a) p- pa\ 'a) 



A6{0,±1} v i,j=l 



(A10) 



(All) 



(A12) 



and 

s = l E { E ^ } + j E (<*« - °») "< • < A13 ) 

A6{0,±1} ^ m,^=0 ^ i=l 

The first term in the r.h.s. of (All) is of hamiltonian form. This is a general 
feature of the reduced dynamics: even in absence of an initial system dynamics, a non- 
trivial hamiltonian contribution is always generated by the dissipative piece (A5). This 
mechanism has deep consequences in neutrino physics: as remarked in the text, it allows 
oscillation phenomena even in case of initially massless neutrinos. 

Further information on the coefficient matrices a^J and b^} in (A9) can be obtained 
by studying the behaviour of the environment correlations in (A10). The operators 
can be taken to be a general linear expression in the environment variables; these are 
the creation, A\(k), and annihilation, A a (k), operators for the quanta representing the 
environment modes, living in an abstract n-dimensional space: 

(A. 14) 

The coefficients xf? "embed" the n-dimensional environment modes into the effective two- 
dimensional neutrino Hilbert space, while f a (k) are appropriate test functions necessary 
to make the operator S M and its correlations well-defined; it can be taken to be of the 
form (|/c|/Mp) m / 2 g a (k), for some positive integer m, with g a (k) of Gaussian form. For 

sake of definiteness, in the following we shall use: g a (k) = g a (n)e- T1 6 /2 , with g a (Q) 
depending only on the angle variables and rj a real constant. The function e(k) gives the 
dispersion relation obeyed by the environment modes; for simplicity we shall adopt an 
ultrarelativistic law: e(k) = \k\ = e. The powers of Mp, characterizing the energy scale of 
the environment, are necessary to give _B M the right dimension of energy. 

We assume an indefinite statistics for these modes, so that the creation and annihila- 
tion operators obey generalized commutation relations: 

A a {k) A\{k') -qAl(k') A a (k) = S ab S^ik - k') ; (A15) 
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the real parameter q determines the mode statistics: the case q = 1 corresponds to stan- 
dard bosons, while for q = one obtains the degenerate algebra discussed in [65-67], in 
connection with DO-branes and black holes. Without loss of generality, we shall assume 
q < 1. Furthermore, the single-mode hamiltonian can be taken to be proportional to the 
corresponding number operator, so that the total environment hamiltonian Hg satisfies 
the relation 



[H £ , Al(k)] = e(k) A\{k) , [H £ , A a (k)} = -e(k) A a (k) 



(A16) 



implicit in the time dependence of (A. 14). 

The thermal correlations involved in the definitions in (A.9) can now be readily com- 
puted. For instance, one explicitly gets: 



„-..> = / J+p-aiAwot / j m+n-3 

2M m+ "" 4 



oo 







where 



a,b,c 



(k) xfxf g c {k) 



2 2 

,-V £ 



(A.17) 
(AA8) 



involves the integration over the angle variables; notice that Xy(e) is a real, symmetric 
matrix. 

This matrix is not generic: it turns out that in order to satisfy the condition of entropy 
increase for finite Xij(e) must vanish for i,j = 1, 2. With this choice, one finds that 

the non- vanishing contributions to L[p] in (A. 12) can come only from the coefficient a 3 ° 3 , 
provided m = 3 — n. Explicitly, one obtains: 



a (0) - 

tt 33 — 



7T 



1 



g 2 M F X 33 (0) . 



(A19) 



The dimensionless coupling constant g, should be expressible in terms of the relevant 
energy scales, i.e. the average neutrino energy E and the mass Mp characteristic of the 
environment. Since g is small, it must be at most of order E/Mp. As a consequence, it 
turns out that the dissipative parameter must scale as E 2 /M F . As mentioned in the 
text, this is a general prediction of the open system approach to dissipation. 

Using the expansion p = as in the text, one immediately finds that the 

dissipative contribution L[p] in (A. 12) is of the form (3.3), with a = a = a 33 and b = c = 
j3 = 7 = 0. In the weak coupling limit a special form of the matrix (3.3) is then selected: 
it is expressible in terms of only one non-standard parameter. This situation does not hold 
any more in the case of the singular coupling limit: in that case, all six parameters a, 6, c, 
a, (3, and 7 are in general non- vanishing (for details, see [21]). 

In a similar way, also the hamiltonian contribution in (A. 13) can be explicitly com- 
puted. Taking into account the results obtained in the evaluation of the coefficients a 
one sees that the 2x2 matrix H becomes diagonal: 

H = a 3 , 



(A) 



(A20) 
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where 



b&> =g 2 M F X 03 (0) G((3 F / V ) , 



(A21) 



and 

G((3 F / V ) = jf dt ^ de sinet " [ 1 _ qe -, F£ ) ■ (^-22) 

In the case of infinite statistics, q = 0, the function G can be explicitly evaluated in terms 
of generalized hypergeometric functions: 

GM^^g; 14) 44)- (A23) 

The operator in (A. 20) contributes via the parameter us = to the effective hamil- 
tonian H e g in (3.1): even in absence of the standard piece u>o, the quantity ujs would 
still generate a level splitting between the two neutrino mass eigenstates, making possible 
oscillation phenomena. 

As a further remark, notice that although both generated via the interaction with 
the environment, the magnitude of the hamiltonian contribution us could differ from the 
dissipative one in (A. 19), since their ratio involves the function G. Although in a different 
context, this phenomenon has also been observed in [57]. 



24 



REFERENCES 



1. R. Alicki and K. Lendi, Quantum Dynamical Semigroups and Applications, Lect. 
Notes Phys. 286, (Springer- Verlag, Berlin, 1987) 

2. V. Gorini, A. Frigerio, M. Verri, A. Kossakowski and E.C.G. Surdarshan, Rep. Math. 
Phys. 13 (1978) 149 

3. H. Spohn, Rev. Mod. Phys. 52 (1980) 569 

4. A. Royer, Phys. Rev. Lett. 77 (1996) 3272 

5. C.W. Gardiner and P. Zoller, Quantum Noise, 2nd. ed. (Springer, Berlin, 2000) 

6. W.H. Louisell, Quantum Statistical Properties of Radiation, (Wiley, New York, 1973) 

7. M.O. Scully and M.S. Zubairy, Quantum Optics (Cambridge University Press, Cam- 
bridge, 1997) 

8. L. Fonda, G.C. Ghirardi and A. Rimini, Rep. Prog. Phys. 41 (1978) 587 

9. G.C. Ghirardi, A. Rimini and T. Weber, Phys. Rev. D 34 (1986) 470 

10. H. Nakazato, M. Namiki and S. Pascazio, Int. J. Mod. Phys. BIO (1996) 247 

11. F. Benatti and R. Floreanini, Phys. Lett. B428 (1998) 149 

12. M.S. Marinov, JETP Lett. 15 (1972) 479; Sov. J. Nucl. Phys. 19 (1974) 173; Nucl. 
Phys. B253 (1985) 609 

13. J. Ellis, J.L. Lopez, N.E. Mavromatos and D.V. Nanopoulos, Phys. Rev. D 53 (1996) 
3846 

14. P. Huet and M.E. Peskin, Nucl. Phys. B434 (1995) 3 

15. F. Benatti and R. Floreanini, Nucl. Phys. B488 (1997) 335 

16. F. Benatti and R. Floreanini, Nucl. Phys. B511 (1998) 550 

17. F. Benatti and R. Floreanini, Phys. Lett. B451 (1999) 422 

18. F. Benatti and R. Floreanini, JHEP 02 (2000) 032 

19. F. Benatti and R. Floreanini, Phys. Rev. D 62 (2000) 125009 

20. J. Ellis, N.E. Mavromatos and D.V. Nanopoulos, Phys. Lett. B293 (1992) 37; Int. 
J. Mod. Phys. All (1996) 1489 

21. F. Benatti and R. Floreanini, Ann. of Phys. 273 (1999) 58 

22. S. Hawking, Comm. Math. Phys. 87 (1983) 395; Phys. Rev. D 37 (1988) 904; Phys. 
Rev. D 53 (1996) 3099; S. Hawking and C. Hunter, Phys. Rev. D 59 (1999) 044025 

23. J. Ellis, J.S. Hagelin, D.V. Nanopoulos and M. Srednicki, Nucl. Phys. B241 (1984) 
381; 

24. S. Coleman, Nucl. Phys. B307 (1988) 867 

25. S.B. Giddings and A. Strominger, Nucl. Phys. B307 (1988) 854 

26. M. Srednicki, Nucl. Phys. B410 (1993) 143 

27. W.G. Unruh and R.M. Wald, Phys. Rev. D 52 (1995) 2176 

28. L.J. Garay, Phys. Rev. Lett. 80 (1998) 2508; Phys. Rev. D 58 (1998) 124015 



25 



29. F. Benatti and R. Floreanini, Phys. Lett. B401 (1997) 337 

30. F. Benatti and R. Floreanini, Phys. Lett. B465 (1999) 260 

31. F. Benatti and R. Floreanini, Irreversibility and dissipation in neutral S-meson decays, 
Nucl. Phys. B, to appear, [hep-ph/0103239| 

32. E. Lisi, A. Marrone and D. Montanino, Phys. Rev. Lett. 85 (2000) 1166 

33. S.M. Bilenky and S.T. Petcov, Rev. Mod. Phys. 59 (1987) 671 

34. C.W Kim and A. Pevsner, Neutrinos in Physics and Astrophysics, (Harwood Aca- 
demic Press, 1993) 

35. M. Fukugita and T. Yanagida, in Physics and Astrophysics of Neutrinos, M. Fukugita 
and A. Suzuki, eds., (Springer, Tokyo, 1994) 

36. R.N. Mohapatra and P.B. Pal, Massive Neutrinos in Physics and Astrophysics, 2nd 
ed., (World Scientific, Singapore, 1999) 

37. S.M. Bilenky, C. Giunti and W. Grimus, Prog. Part. Nucl. Phys. 43 (1999) 1 

38. J. Schechter and J.W.F. Valle, Phys Rev. D 22 (1980) 2227; ibid, bf 23 (1981) 1666 

39. L.F. Li and F. Wilczek, Phys Rev. D 25 (1982) 143 

40. B. Kayser, S.M. Bilenky, J. Hosek, S.T. Petcov, Phys. Lett. B94 (1989) 495; M. Doi, 
T. Kotani, H. Nishiura, K. Okuda and E. Takasugi, Phys. Lett. B102 (1981) 323; P. 
Langacker et al. Nucl. Phys. B282 (1987) 589 

41. F. Benatti and R. Floreanini, Mod. Phys. Lett. A12 (1997) 1465; Banach Center 
Publications, 43 (1998) 71; Phys. Lett. B468 (1999) 287; On the weak-coupling limit 
and complete positivity, Chaos Sol. Frac, to appear 

42. T.K. Kuo and J. Pantaleone, Rev. Mod. Phys. 61 (1989) 937 

43. S.P Mikheyev and A.Yu Smirnov, Sov. J. Nucl. Phys. 42 (1985) 913; Nuovo Cim. 
9C (1986) 17 

44. L. Wolfenstein, Phys. Rev. D 17 (1978) 2369; ibid. 20 (1979) 2634 

45. For instance, see: M. Artin, Algebra, (Prentice Hall, Englewood Cliffs (NJ), 1991) 

46. K. Lendi, J. Phys. A 20 (1987) 13 

47. A.M. Gago, E.M. Santos, W.J.C. Teves, R. Zukanovich Funchal, Quantum dissipative 
effects and neutrinos: current constraints and future perspectives, hep-ph/ 0009222 

48. S. Adler, Phys. Rev. D 62 (2000) 117901 

49. J. Ellis, N.E. Mavromatos and D.V. Nanopoulos, Phys. Rev. D 63 (2001) 024024 

50. Y. Liu, L. Hu and M.-L. Ge, Phys. Rev. D 56 (1997) 6648 

51. C.-H. Chang, W.-S. Dai, X.-Q. Li, Y. Liu, F.-C. Ma and Z. Tao, Phys. Rev. D 60 
(1999) 033006 



52. CP. Sun and D.L. Zhou, Quantum decoherence effect and neutrino oscillations, [hep- 
ph/980833^ 



53. A.D. Dolgov, Sov. J. Nucl. Phys. 33 (1981) 700 

26 



54. G. Raffelt, G. Sigl and L. Stodolsky, Phys. Rev. Lett. 70 (1993) 2363; G. Sigl 
and G. Raffelt, Nucl. Phys. B406 (1993) 423; G. Raffelt, Stars as Laboratories for 
Fundamental Physics, (Chicago University Press, Chicago, 1996) 

55. M. Prakash, J.M. Lattimer, R.F. Sawyer and R.R. Volkas, Neutrinos in dense matter, 
early universe, supernovae, neutron stars, |atro-ph/ 0103095 

56. T. Ohlsson, Equivalence between Gaussian averaged neutrino oscillations and neutrino 
decoherence, [tiep-ph/00 12272 

57. CP. Burgess and D. Michaud, Ann. of Phys. 256 (1997) 1; P. Bamert, CP. Burgess 
and D. Michaud, Nucl. Phys. B513 (1998) 319 

58. P. Lipari and M. Lusignoli, Phys. Rev D 60 (1999) 013003 

59. G.L. Fogli, E. Lisi, A. Marrone and G. Scioscia, Phys. Rev. D 60 (1999) 053006 

60. J. Alfaro, H.A. Morales- Tecotl and L.F. Urrutia, Phys. Rev. Lett. 84 (2000) 2318 

61. C.N. Leung, Neutrino tests of general and special relativity, |hep-ph/0002073| 

62. H.V. Klapdor-Kleingrothaus, H. Pas and U. Sarkar, Effects of quantum space time 
foam in the neutrino sector, |hep-ph/ 0004123] 

63. V. Barger, S. Pakvasa, T.J. Weiler and K. Whisnant, CPT-odd resonances in neutrino 
oscillations, piep-ph/ 0005197 



64. M. Lindner, T. Ohlsson and W. Winter, A combined treatment of neutrino decay and 
neutrino oscillations, [hep-ph/0 103170] 

65. A. Strominger, Phys. Rev. Lett. 71 (1993) 3397 

66. I. V. Volovich, D-branes, black holes and SU (oo) gauge theory, |tiep-th/9608137 



67. D. Minic, Infinite statistics and black holes in Matrix theory, Pennsylvania University 
preprint, 1997, |hep-th/9712202 



27 



